A model independent determination of the 

B — > X s j decay rate 



Proceedings of CKM 2012, the 7th International Workshop on the CKM Unitarity 
Triangle, University of Cincinnati, USA, 28 September - 2 October 2012 



Florian U. Bernlochner 

University of Victoria, Victoria, British Columbia, Canada V8W 3P 
E-Mail: florian@slac. Stanford, edu 

Heiko Lacker 

Humboldt University of Berlin, 12489 Berlin, Germany 
E-Mail: lacker@physik.hu-berlin. de 

Zoltan Ligeti 

Ernest Orlando Lawrence Berkeley National Laboratory, University of California, 
Berkeley, CA 94720, USA 
E-Mail: ligeti@lbl.gov 

Iain W. Stewart 

Center for Theoretical Physics, Massachusetts Institute of Technology, 
Cambridge, MA 02139, USA 
E-Mail: iains@mit.edu 

Frank J. Tackmann 

Deutsches Elektronen- Synchrotron (DESY), D-22607 Hamburg, Germany 
E-Mail: frank.tackmann@desy.de 

Kerstin Tackmann 

Deutsches Elektronen- Synchrotron (DESY), D-22607 Hamburg, Germany 
E-Mail: kerstin. tackmann© desy. de 

MIT-CTP 4429, DESY 13-019 



1 Introduction 

Rare penguin decays play an important role in the search for new physics in the 
flavor sector: the presence of new heavy particles which couple to heavy quarks 




Figure 1: On the left-hand side the SM penguin B — > X s "y decay is shown. In 
new-physics models, such as e.g. the minimal supersymmetric Standard Model, the 
Wilson coefficient associated with the effective operator Or shown at the right-hand 
side is modified. 

would manifest itself in modifications of the total decay rate. In this presentation, a 
model-independent analysis is carried out to test the compatibility of four B — > X s 7 
measurements with the Standard Model expectation (SM). In Fig. [TJthe SM decay 
and a possible extension via a charged Higgs boson are shown. In the past, such 
an analysis used the extrapolation of the measured decay rates to the low E 1 region 
to compare the experimental measured partial branching fraction with the next-to- 
next-to-leading order SM prediction, B{E 1 > 1.6 GeV) = (3.55 ± 0.24 ± 0.09) x 10" 4 , 
from Ref. [U |2], as adopted e.g. by the analysis of Ref. [3] and more recent updates. 
The reason for this extrapolation lies in the theoretically poorly known nonpertur- 
bative corrections from the b quark distribution function in the B meson, called the 
shape function, which affects the partial decay rate, B(E^ > E^ ut ), for high values of 
E^ ut . The drawback of the extrapolation based analyses lies in the introduction of 
undesired model dependencies, which are hard to quantify, and also in making poor 
use of the experimentally most precise regions at high E 1 by including the partial 
branching fraction down to the region of phase space with low E 1 , which is domi- 
nated by large background contributions from other B meson decays. The SIMBA 
collaboration uses an alternative approach, outlined in Ref. by determining the 
shape function directly from data from a global analysis of the B — > X s 7 spectra with 
negligible model dependence compared to the present experimental and theoretical 
uncertainties. The main changes with respect to our earlier work Ref. [5j |6] lies in 
the evaluation of the uncertainties of missing higher-order perturbative corrections, 
and the inclusion of the measured E 1 spectrum of Ref. [7j. 
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2 Treatment of the shape function and formulae 
for B — > X s 7 decay rate 

2.1 Shape function 

The shape function, which enters in the differential B — > X s 7 decay rate, can be 
factorized into perturbative and non-perturbative contributions, cf. [4], via 



S(co, h) = J dk C {co - k, fi) F{k) 



where Cq{u) — k, jj) is the MS-renormalized 6-quark matrix element of the shape func- 
tion operator calculated in perturbation theory, and F(k) are the non-perturbative 
contributions to S(u,(jl). Constructing the shape function as done in Eq. [I] offers 
several advantages over alternative approaches: the shape function has the correct 
perturbative tail at large u, and the correct RGE behavior; for small u, the shape 
function is dominated by the non-perturbative parameter F(k). This means that the 
shape of the B — > X s 7 spectrum at large E 1 is determined by the non-perturbative 
parameter F(k). 

The approach outlined in Ref. |4] proposed the determination of F(k) directly from 
experimental spectra by employing an expansion in a set of complete orthonormal 
basis functions, f n , as 



A 



Cn f n (a) 



n=0 



/oo 
dkF(k) = J2°n =!> ( 2 ) 
n=0 



where A ~ Aqcd is a dimensional parameter of the basis. Since the orthonormal 
basis in Eq. [2] is complete, this description offers a model-independent description 
of the shape function and the B — > X s ^ decay rate, and the shape of the spec- 
trum is parametrized by the expansion coefficients c n . These coefficients can be 
determined directly by fitting the available data, taking into account the full exper- 
imental uncertainties and correlations. In practice, however, the expansion in Eq. [2] 
needs to be truncated, since the available experimental information only allows for 
a finite number of N coefficients to be constrained. This truncation introduces a 
residual model-dependence, which depends on the chosen functional basis and scales 
as 1 — J2n=o c n- The optimal values for A and N need to be determined from data: A 
is chosen such that the fitted series converges quickly, and the number of used basis 
functions N should be chosen large enough such that the truncation uncertainty is 
small compared to the experimental uncertainties. 
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2.2 Master formulae for B — > X s ^f decay rate 
The B — > Xg'-f photon energy spectrum is given by 

dr G 2 F a em ^3^2 it/ T /*|2 
E m b \ VtbV ts \ 



d£ 7 2tt 4 



x ||Crf / d& W 77 (k)F(m B - 2£ 7 - fc) + J] W 77>m F m (m B - 2£ 7 ) 

dA; ^ [2Re(C 7 ncl )a W^fc) + C i C j %{k)^ F(m B - 2E 1 - k) j , (3) 



+ 



where H^7 7 (/c) contains the perturbative corrections to the b — > 57 decay via the elec- 
tromagnetic dipole operator, O7, resummed to next-to-next-to-leading-logarithmic or- 
der [Hill], and including the full NNLO corrections [91 lioj^j F[k) is the leading shape 
function as introduced in the previous section, the F m {k) denote l/fhb suppressed 
subleading shape functions, and W 7i (k) and Wij(k) contain subleading perturbative 
corrections. The full set of expressions entering Eq. [3] will be given in Ref. [TT] . In 
a fit to B — » X s 7 spectra the subleading shape functions can be absorbed at lowest 
order in a s into the leading shape function, reducing the number of coefficients which 
need to be determined from data. The coefficient C 7 ncl multiplying the dominant 
contribution proportional to W 77 (k) in Eq. |3l is defined as 



Of* = Cf + £ C^ ) + r 8 M Cf + . . . , (4) 

i=i 

where Cf s (fio) are the standard scheme-independent effective Wilson coefficients and 
fnb(no) is the MS 6-quark mass. The coefficients ri_6,s(A t o) contain all virtual contri- 
butions from the operators Oi_6,s that generate the same effective b — > 37 vertex as 
O7. The ellipses denote included terms proportional to ln(/io/m b ) that are required 
to cancel the /x dependence on the right-hand side and vanish at /i = fh b , such that 
C 7 ncl is by definition /io-independent to the order one is working at. 

Since the terms in the last line in Eq. [3] are small, we consider |C 7 ncl VtbV t * | as 
the parameter that determines the normalization of the B — > Xgj rate. This nor- 
malization is extracted simultaneously with F(k) from a fit to the various measured 
Ej spectra. The important contributions from Oi_6,s are the virtual corrections con- 
tained in C 7 ncl , which have a sizable effect on the normalization of the B — > X s ^f rate. 
By including them in C 7 ncl , they explicitly do not affect the shape of the spectrum, 
and so do not enter in our fit. They instead enter in the SM prediction for C 7 ncl , which 
can be computed independently. Below, the fit result is compared to the NLO SM 

1 At lowest order in perturbation theory: Wn(k) = 8{k) 
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value, Cf cl = 0.354±g;gi| [II]. For a more stringent test for new physics, evaluating 
C^ cl in the SM at NNLO along the lines of Refs. [U [2] would prove very valuable. 

3 Fit to the available B-Factory data 

3.1 Experimental data from BaBar and Belle 

As experimental inputs the Belle measurement from Ref. [12], and the three BaBar 
measurements from Refs. [131 [7J are used. The experimental statistical and sys- 
tematic uncertainties and correlations are fully included in the fit procedure. The 
BaBar spectra of Refs. [T3l E] are measured in the B rest frame and are corrected 
for efficiencies. The experimental resolution in E 1 for each spectrum is smaller than 
its respective bin size, so both spectra can be directly used in the fit. The Belle spec- 
trum from Ref. [12] is measured in the T(4S) frame and affected by both efficiency 
and resolution. The BaBar result of Ref. [7J provides resolution unfolded spectra 
with smaller bin- by-bin correlations in the B and the T(4S f ) frame. We analyze the 
unfolded result in the T(45*) frame. All four spectra are shown in Fig. [2j 

3.2 Fit setup 

To fit to the experimentally measured photon energy spectra, the expansion for F(k) 
in Eq. [2] is used for Eq. [3] and integrated over the appropriate range of E 1 for each 
experimental bin and each combination of basis functions f m {x)f n (x). The theory 
prediction for the ith bin, B l , is then given by 

N 

B* = ml \Cr l V tb V;f CmC n B* mn + ■■■, (5) 

m,n=0 

where the ellipses denote the additional included terms arising from the last line in 
Eq. [3J The overall ttt| is expressed in terms of the moments of F(k), so it is effectively 
a function of the c n coefficients. A \ 2 minimization is performed to all available bins 

with cox-,n and \C^ cl VtbV t * s | as the fit parameters. The constraint CqH Vc 2 N = 1 is 

enforced to ensure that F(k) is properly normalized to unity. This fitting procedure 
was extensively tested using pseudo-experiments and provides unbiased central values 
with correct uncertainties. 

3.3 Fit results 

For our default fit a value of A = 0.5 GeV and four basis coefficients 00,1,2,3 was chosen: 
These values were selected after checking carefully the convergence of various basis 
expansions, number of coefficients, and with the general focus of avoiding over-tuning. 
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Figure 2: Used Belle and BaBar B — > X s ^y measurements: Ref. [13] (top left), 
Refs. [7j (top right), Ref. [H] (bottom left), Ref. [12] (bottom right); The histograms 
show the result of the fit with a basis of A = 0.5 GeV with two (yellow), three (green), 
four (blue), and four (orange) coefficients. The default fit result uses four coefficients. 

The corresponding fits with two (co,i), three (00,1,2), four (00,1,2,3), and five (00,1,2,3,4) 
expansion coefficients with A = 0.5 GeV are shown in Fig. [2] The fit converges 
after the inclusion of four coefficients and describes the measured spectra well. The 
X 2 /ndf for the default fit with four coefficients is 41.65/48 corresponding to a p- value 
of 0.87. The fit results for the shape function for 2,3,4, and 5 basis coefficients are 
shown in the left panel of Fig. |3j The corresponding results for | Cy ncl Vf bV^.* | and 
m\ s , where the latter is computed from the moments of the fitted F(k), are shown 
in the right panel of Fig. [3j The shape function in Fig. [3] verify the convergence of 
the basis expansion as the number of basis functions is increased. As one expects, 
the uncertainties returned by the fit increase with more coefficients due to the larger 
number of degrees of freedom. However, with too few coefficients one would have to 
add the truncation uncertainty. A reliable value for the final uncertainty is provided 
by the fitted uncertainty when the central values have converged and the respective 
last coefficients, here C3 or C4, are compatible with zero. At this point, the truncation 
uncertainty can be neglected compared to the fit uncertainties. Equivalently, the 
increase in the fit uncertainties from including the last coefficient that is compatible 
with zero effectively takes into account the truncation uncertainty. Using a fixed 
model function and fitting one or two model parameters would thus underestimate 
the true model uncertainties in the shape function model. 



6 



3.4 Theoretical uncertainties from missing higher order per- 
turbative corrections 

Including the theory uncertainties from missing higher order perturbation corrections 
for the default fit gives 

\Cr l V tb V*\ = (14.83 ± 0.53 [cxp] ± 0.37 [thco] ) x 1(T 3 , 

m\ s = 4.77 ± 0.03 [cxp] ± 0.02 [theo] , (6) 

with comparable sizes for experimental and theoretical uncertainties. The result for 
|C} ncl | is compatible within one sigma with the NLO SM value, shown as a grey band 
in Fig. [3} for which a value of |Vtf>14*| = 40.68±°;5 x 10 -3 was used. 

4 Summary and outlook 

We presented preliminary results from a global fit to B — >■ X s 7 data, which deter- 
mines the total B — > X s ^ rate, parametrized by |C7 ncl VtfcVJ*|, and the S-meson shape 
function within a model-independent framework. The value of |Cy ncl VtfeV^* | extracted 
from data agrees with the SM prediction within uncertainties. From the moments 
of the extracted shape function we determine m\ . In the future, information on 
rrib from other independent determinations can be included. The shape function 
extracted from B — > Xg'j is an essential input to the determination of \V u b\ from 
inclusive B — > X u tv decays. 

A combined fit to B — > Xg'j and B — > X u lv measurements within our framework 
is in progress. It will allow for a simultaneous determination of \Cy V^V^*] and 
\V u b\ along with the shape function with reliable uncertainties. In addition to a 
few branching fractions with fixed cuts, it is important to have measurements of 
the B — > X u tv differential spectra (including correlations), e.g. the lepton energy 
or hadronic invariant mass spectra. As for B — > X^, fitting the differential spectra 
allows making maximal use of the B — > X u lv measurement, by using the experimental 
most precise regions to constrain the nonperturbative inputs and further reduce the 
associated uncertainties. 
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Figure 3: The extracted F{k) (with absorbed l/m b corrections) for two (c 0) i), three 
( c o,i,2)) f° ur ( c o,i,2,3)) ano - five (00,1,2,3,4) coefficients and basis parameter A = 0.5 GeV 
(Left). The colored envelopes are given by the uncertainties and correlations of the 
extracted coefficients c n . Determined values of \C^ cl V t bV t *\ and m\ s with A% 2 = 1 
contour (Right). The grey band shows the SM value of |C*7 ncl V^f,V^*| and the pink 
band show the SCET scale variation uncertainties. 
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